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ABSTRACT 

We present an occultation of the central X-ray emitting region in the Seyfert 
Galaxy NGC 1365. This extreme spectral variation (from Compton-thin to 
reflection-dominated and back to Compton-thin in four days) has been caught in 
a ten days Chandra monitoring campaign consisting of six short (15 ks) observa- 
tions performed every two days. We discuss the implications of this occultation 
within the scenario of a Compton-thick cloud crossing the line of sight of the 
X-ray source. We estimate a source size R < 10 14 cm and a distance of the 
cloud from the source D < 10 16 cm. This direct measurement confirms the the- 
oretical expectations of an extremely compact X-ray source, and shows that the 
Compton-thick circumnuclear gas is located at a distance from the center on the 
scale of the Broad Line Region. 

Subject headings: Galaxies: AGN — Galaxies: individual (NGC 1365) 

1. Introduction 

Variability of the X-ray absorbing gas is common in Active Galactic Nuclei on time 
scales from months to years (Risaliti, Elvis & Nicastro 2002). If an observed variation of 
absorbing column density Nh is due to clouds in virialized motion crossing the line of sight, 
the amount and the duration of the Nh variation put constraints on the distance of the 
obscuring cloud from the center, and on the density of the cloud. 

In a few cases an Nh variation has been detected within a single observation (NGC 4388 
Elvis et al. 2004, NGC 4151, Puccetti et al. 2007), indicating that the absorber must be 
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extremely compact, i.e. on the scale of or slightly larger than the Broad Line Region (~ 
10 16 cm for a 10 8 M black hole). 

In this context, the Seyfert Galaxy NGC 1365 has shown extreme variability in the past 
~ 12 years: it was observed in a reflection-dominated state by ASCA in 1995 (Iyomoto et 
al. 1997), then in a Compton-thin state by BeppoSAX in 1997, with absorbing column density 
Nh ~ 4 x 10 23 cm -2 (Risaliti et al. 2000). Such a long time interval between the observations 
leaves two possible scenarios open: extreme absorption variability (A(iV#) > 10 24 cm~ 2 ) or 
a switch off and on of the X-ray source. This ambiguity has been solved with a more recent 
set of short observations, performed by Chandra in December 2002, and by XMM-Newton 
three and six weeks later. The source was caught in a reflection-dominated state in the 
first and third observation, while it was in a Compton-thin state with N H ~ 4 x 10 23 cm -2 
in the second observation (Risaliti et al. 2005A, hereafter R05). Such fast variations are 
hard to explain within the intrinsic variation scenario strongly suggesting that the observed 
reflection-dominated states are due to Compton-thick clouds crossing the line of sight (see 
R05 for a full discussion). Three additional XMM-Newton observations performed in 2003 
and 2004 caught the source in a Compton-thin state, with column densities between 1.5 and 
5xl0 23 cm" 2 . 

The latest two observations are relatively long (60 ks) and allowed a detailed timing and 
spectral analysis. The spectral analysis revealed the presence of a highly ionized, compact 
absorber (Risaliti et al. 2005B), while the timing analysis revealed column density variability 
of A(N H ) ~ 10 23 cm" 2 on time scales of ~ 50 ks (Risaliti 2006, Risaliti et al. 2007, in prep.). 

In order to explore the variability time scales between the longest single observations 
(~a day) and the shortest observed Compton-thick/Compton-thin change (three weeks) we 
conducted a Chandra campaign consisting of six 15 ks observations performed once every 
~ 2 days for ten days in April 2006. Here we report the results of these observations, with 
emphasis on an occultation event which occurred during the first four days of monitoring, 
and we discuss the physical implications of these results. 

2. Reduction and Data Analysis 

The observation log is shown in Tab. 1. All the observations were performed with the 
ACIS-S instrument (Weisskopf et al. 2002) in T/4 window' mode in order to avoid possible 
pile-up. A check of the reduced spectra confirmed that in all cases the pile-up is lower than 
1%. 
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The data were reduced using the CIAO 3.3(11 package and using a standard procedure, 
as described in the CIAO Threads 1 . We extracted the spectrum from a circular region with 
a 2 arcsec radius. This removes most of the soft, diffuse emission from the spectra (see Fig 1 
in R05). A complete analysis of the spectral and spatial properties of the diffuse component 
will be presented elsewhere (Bianchi et al. 2007, in prep.). The background was selected from 
a region in the field free from contaminating sources. The spectral analysis was performed 
using the SHERPA package inside CIAO. 

A simple visual inspection of the spectra (Fig. 1) is enough to notice the main result 
of this work: the second spectrum is completely different from the others both in flux (it is 
much fainter) and in shape (it is flatter, with a prominent emission feature at ~ 6.5 keV), 
indicating that the source switched to a reflection-dominated state between the first and 
second observation, and then switched back to a transmission-dominated state between the 
second and the third observation. 

We first performed a simple spectral analysis, with a model consisting of a thermal 
emission at low energies, an absorbed power law and an iron emission line with E=6.4keV. In 
all the fits discussed below, the soft component is fitted with a thermal component (Raymond 
& Smith 1977) with temperature kT=0.8 keV and constant within 2% in all cases. 

All the high energy (E=2-10 keV) spectra except for the second one are well fitted 
(reduced \ 2 between 1.1 and 1.5) with a power law with photon index T = 1.8 — 2, a column 
density Njj = 2 — 4 x 10 23 cm -2 and an iron line with equivalent width EW=200-300 eV. 
The second spectrum is fitted by a flat power law with T ~ 0.5 and Njj < 10 21 cm -2 . The 
emission feature has an equivalent width EW~ 1.5 keV. 



1 http:/ / asc.harvard.edu/ciao/ 
Table 1: NGC 1365 - New Chandra Observations 
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Observation log of the six Chandra ACIS-S observations, started on April 14, 2006. a : Time 
elapsed from the end of the previous observation. 
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These results confirmed the visual inspection, and prompted a more detailed analysis, 
which has been performed in two main steps: 

1) We fitted the second spectrum with a cold reflection continuum (PEXRAV model, Magzdi- 
arz & Zdziarski 1995) and an emission line at 6.4 keV, corresponding to neutral iron Ka emis- 
sion. A second line at energy E=6.9 keV, corresponding to hydrogen-like iron, is also required 
by the fit. The photon index of the intrinsic component is not well constrained, so we fixed 
it to the average of the values obtained fitting the transmission-dominated spectra. The best 
fit equivalent widths of the two lines are EW 6 . 4 = 1.2 ± 0.2 keV and EW 6 . 9 = 0.9 ± 0.4 keV. 

2) The other five spectra were fitted with the same reflection continuum (with all the pa- 
rameters frozen to the best fit values obtained in the analysis of the second spectrum) plus a 
continuum component representing the direct emission of the X-ray source. Since the spectra 
obtained in the past XMM-Newton observations have a much higher signal-to- noise, we used 
the best fit models obtained from those observations as a baseline for our spectral fitting. 
The best fit model consists of an absorbed power law, plus four narrow absorption lines 
between 6.7 and 8.3 keV (representing He-like and H-like iron Fe Ka and Kf3 transitions), 
and a broad emission line. We refer to Risaliti et al. 2005B for a detailed description of these 
models. 

The results of the fits are in all cases satisfactory from a statistical point of view (reduced 
X 2 between 1 and 1.2) and provide best-fit values for the continuum photon index constant 
within the errors, while the values of the column density show significant variations. We 
then repeated the whole analysis fitting all the spectra simultaneously, requesting a constant 
value for the photon index and leaving the other parameters free to vary. We obtained as 
good a fit as in the previous case from a statistical point of view (overall reduced x 2= l-01) 
and slightly smaller error intervals for the column density estimates. As a side product of our 
analysis, we mention that the two strongest iron absorption lines (K a lines from He-like and 
H-like iron) are significantly detected. The inclusion of these features in the fit is relevant 
in our context only because neglecting them would affect the estimate of the continuum 
parameters. However, their detection is an interesting independent confirmation of the highly 
ionized absorber discovered in the XMM-Newton observations (Risaliti et al. 2005B). We will 
discuss this issue elsewhere (Risaliti et al. 2007, in prep.). 

The main best fit parameters are shown in Table 2. 

3. Discussion 

The analysis described in the previous Section shows that a change from Compton-thin 
to Compton-thick states, and then back to Compton-thin, occurred in the first four days of 
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Fig. 1. — Spectra obtained from the first three Chandra observations of NGC 1365. the 
time interval between the single observations is two days. The first and third spectrum are 
typical of a transmission-dominated source, with a steep continuum and a photoectric cut-off 
at ~ 4 keV. The second spectrum is fainter by a factor ~ 10, and is characterized by a flat 
continuum, and a prominent iron emission line, typical of reflection-dominated sources. The 
last three Chandra observations are similar to the first and third ones, and are not shown 
for clarity. 
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Fig. 2. — Schematic representation of the observed eclipse. The intervening thick cloud 
(orange circle, diamater Dc) vith Keplerian velocity Vk starts covering the X-ray source 
(blue circle, diameter D$) at some time between the first and second Chandra observations. 
After a time T\— Ds/Vk the source is completely covered, and remains obscured for a time 
T2 = {Dc — Ds)/Vk- In this state it is observed for the second time by Chandra. Then, it 
gradually uncovers, until it is back in the initial state. After some more time it is observed 
again by Chandra. From the times between the observations we infer T 2 + 27\ < 4 days; 
T 2 > 4 hours. The lower part of the figure shows the time evolution of the observed flux. 
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our Chandra monitoring of NGC 1365. 

As discussed in R05, such rapid changes are hard to explain with intrinsic spectral 
variations. In order to further check this possibility, we estimated the upper limit of a 
possible direct component in the reflection-dominated spectrum. Assuming a spectral shape 
analogous to that obtained for the first and third observations (which are quite similar both 
in flux and spectral parameters, Table 2) the 90% upper limit to the direct flux is only 5% 
of that of observations one and three. A decrease of the intrinsic flux by a factor of 20, and 
then an increase back to the initial flux, besides implying an unlikely "fine tuning" in order 
to reproduce the observed symmetry between the fading and recovery phase, would require 
a cooling time of at least 3-4 weeks in the framework of a Shakura-Sunyaev (1973) disk (see 
R05 for details), completely incompatible with our observed varation times. For this reason, 
in the following we will only discuss the occultation scenario. 

In the simplest scheme, shown in Fig. 2, the size of the X-ray emitting source D$ is 
given by the obscuring cloud velocity, Vk, times the ingress/egress time T\. The linear size 
of the obscuring cloud is D c > D s and the distance between the cloud and the source is 
R. Since our goal is to put an upper limit on the source size, it is particularly important to 
discuss the upper limits on the estimates of these two parameters. 

Here we discuss several constraints on Vk, T\ and D S - 
A. Statistical limits on 7\. The occultation event is characterized by two times (Fig. 2): 
the ingress/egress time 7\ = D s /Vk and the time during which the source remains com- 
pletely obscured, T 2 = (D c — D s )/Vk- Considering the X-ray observational history of 
NGC 1365, we note that the source has been observed four times in a completely reflection- 
dominated state, while it has never been observed during the ingress/egress phase, in any 
published observation. In order to check this, we re-analyzed all the past observations look- 
ing for (a) fast drops in the light curves of the transmission-dominated spectra, indicative of 
a possible occultation event during the observation, and (b) possible direct continuum com- 
ponents in the reflection-dominated spectra. The first check easily ruled out the possibility 
that such occultations happened during any observation (this would imply a drop by a factor 
of at least 10 in the hard X- ray (E>4 keV) light curve, which would be easily detected). 
In the second check we added a direct continuum component to the model of the reflection- 
dominated spectra, consisting of an absorbed power law. We required that the photon 
index and the absorbing Nh varied within the lowest and highest values measured in all the 
transmission-dominated spectra. In the three most recent reflection-dominated spectra this 
extra component is not required, and the upper limit to the flux of this component is as low as 
10% of the faintest transmission-dominated spectrum. The result for the ASCA observation 
is inconclusive, because of the lower signal-to- noise. Also in this case the extra component is 



- 8- 



not statistically required, but a significant contribution cannot be ruled out. Therefore, we 
do not include this observation in our analysis. Summarizing, in three cases the source was 
found in a completely reflection dominated state, and in no cases it was caught during an 
ingress/egress. During an eclipse the source is in the partial occultation phase (state 1) for 
a total time 2Ti, and is totally covered (state 2) for a time T 2 . The fraction of time in the 
completely obscured state is / = T 2 /Ttot where Ttot = (2Ti + T 2 ) < 4.2 days based on our 
observations. We require that the probability of finding the source three times in state 2 but 
never in state 1 is P=/ 3 >10%. This implies T 2 > 0A6T TOT and T x < 0.27 x T TOT < 1 day. 
Furthermore, we get T 2 > 1.771, and D c > 2.7D S . 

B. Physical limits on V K . A first limit on V K comes from the measured width of the 
iron emission line in the reflection-dominated spectra. Here we assume that the reflecting 
Compton-thick gas is the same one which is responsible for the occultations. This is in agree- 
ment with the available statistics: since during the whole observational history of the source 
the occurrence of Compton-thick states is 4/14, it is expected that the obscuring clouds 
cover a significant fraction, ~l/3, of the solid angle as seen from the source, and therefore 
they are also expected to contribute significantly to the observed reflection (a typical fitting 
value of the PEXRAV parameter R, normalized to 1 for a half-solid angle coverage of the 
reflector, is indeed ~ 0.6 — 0.8, as shown in Tab. 2). The best available estimate of the iron 
line width, Wfc comes from the XMM-Newton observation of the reflection-dominated state 
of NGC 1365. We obtain W Fe < 150 eV, corresponding to a velocity V <7,000 km/s. The 
line width measures only the average line-of-sight velocity, while during the occultation the 
cloud velocity vector lies in the plane of the sky. Assuming circular orbits, this implies that 
the actual transverse velocity during occultation can be as large as (n/2) x V =12,000 km/s. 
Finally, using T\ < 1 day we obtain D s < 10 14 cm. 

C. Geometrical limits on D s . A geometrical limit on D s can be obtained assuming 
the minimum possible distance R between the source and the cloud. This is given by 
R = (D s + D c )/2. Since the cloud size D c must be equal to or larger than the source size D s , 
the minimum distance for a given Ds is R m in = Ds- If the cloud is moving with Keplerian 
velocity V K = (GM BH /D S ) 1/2 the condition V K x T\ = D s implies D s = (GM BH y/ 3 T 2/3 . 

Two independent estimates of the black hole mass in NGC 1365 are available: \og(M BH /M & ) = 
7.3 ±0.4(0.3) from the M^-bulge velocity dispersion correlation (Oliva et al. 1995, Ferrarese 
et al. 2006) and \og{M BH /M & ) = 7.8 ± 0.4(0.3) from the relation between M BH and the K 
magnitude of the host bulge (Dong & De Robertis 2006, Marconi & Hunt 2003) where the 
errors include statistical and systematic effects, and the number in brackets refer to the 
statistical dispersion of the correlation. 

Using Ti < 1 day we obtain D s <3x 10 14 (M BH /M BES t) 1/3 cm, where M BEST is the 
average estimate of the black hole mass according to the correlations cited above. Expressing 
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the source size in units of gravitational radii we find D s < 33(M BH / 'M B est)~ 2 ^ Rg- We 
note that the higher the black hole mass, the larger the physical size of the source, and 
the smaller the source size in units of gravitational radii. Considering the relatively high 
uncertainty in the black hole mass determination, we can repeat the above calculations using 
the minimum and maximum values compatible with the two correlations, logM B H( m in) = 7. 1 
and \ogM BH ( max } = 8.0. We obtain Dsijnin) < 2 x 10 14 cm, corresponding to 56 Rg, and 
D s (max) < 4.6 x 10 14 cm, corresponding to 15 R G . 

D. Physical limits on R. A conceptually different limit on the distance R between the 
source and the obscuring cloud can be obtained by considering the ionization state of the 
latter. This is defined by the ionization parameter, U = L x /(nR 2 ), where n = N H /D C) is 
the density of the obscuring cloud, located at a distance R from the central source, whose 
X-ray total luminosity is L X - A limit on the ionization parameter can be obtained from the 
analysis of the reflection-dominated spectra observed with XMM- Newton and Chandra: we 
added to the best fit model an extra component consisting of the continuum observed in 
transmission-dominated states, absorbed by gas with variable column density and ionization 
parameter (Done et al. 1992). Since the upper limit on the flux of the direct emission is only a 
few percent of that observed in transmission-dominated states, the absorber must be effective 
enough to remove it almost completely. This implies a lower limit on the absorbing column 
density N B > 10 24 cm -2 and an upper limit on the ionization parameter U < U max = 100. 
From the latter limit, assuming that the cloud is moving with Keplerian velocity Vk and 
requiring that the cloud dimension Dq is larger than the source dimension Dg we obtain, after 
a little algebra: R > (GM B ^) 1//5 [(T 1 + T 2 )L x /(U max N H )} 2 ^. Assuming fiducial values for 
the occultation times and the cloud column density (Ti + T 2 ~ 2 days and N B = 10 24 cm -2 ) 
we obtain R > 3 x 10 15 \M BH /M BE st) 2 ^ cm. We note that the dependence on the black 
hole mass is weak, so even adopting the extreme values allowed for M BB would change the 
result by a factor of less than 30%. From this estimate we can easily obtain a new limit on 
the Keplerian velocity V K = (GM BH /R) 1 ' 2 < 12,000(M BH /M BEST ) 2 / 5 km/s, i.e. the same 
upper limit as using the previous argument on the line width. Analogously, the limit on the 
source size is D s = V K Ti < 10 14 (M BH / mbest) 2/5 cm. 

Summarizing the constraints on the source size, we obtained: 

1) D s < 3 x 10 14 (M S tf/M B £5 T ) 1//3 cm (from geometrical considerations) 

2) Ds < 10 14 cm (from the observed limits on Vk) 

3) D s < 10 14 (M BH I M BE st) 2 ^ cm (from limits on the ionization state of the obscuring 
cloud) 

All these independent methods strongly suggest that the X-ray source has a size of the order 
of or smaller than 10 14 cm, corresponding to 10 R G for a black hole mass M BH = M BEST = 
3 x 10 7 M m . 
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We note that our estimates can be slighly altered if the case of a partial covering of 
the thick cloud in the first and third observations is considered. Specifically, in order to 
reproduce the observed symmetry in the first three observations, the same fraction F of the 
source should be covered in observations one and three. However, if the fraction F is small, 
the effect is negligible, while if F is big, the cloud would uncover the source in the subsequent 
observations, implying a large increase of the flux from the third to the fourth observations, 
which is not observed (Table 2). We conclude that the possible effect of partial covering of 
the transmission-dominated spectra would not significantly affect our estimates of the source 
and cloud sizes. 

Using the calculations discussed above we can easily estimate the distance of the ob- 
scuring cloud from the center. Assuming V K =12,000 km/s and M BEST = 3 x 1O 7 M we 
obtain R = 3 x 10 15 cm, compatible with the lower limit deduced from the ionization con- 
dition. Such a distance corresponds to 300 Rg, and is of the same order of the distance 
from the center of the innermost broad line clouds, and therefore much more compact than 
the circumnuclear medium assumed in the standard Unified Models of AGNs (e.g. Krolik & 
Begelman 1988). 

This work has been partially funded by NASA Grant G06-7102X. We also acknowledge 
financial contribution from contract ASI-INAF 1/023/05/0. 

REFERENCES 

Done, C, Mulchaey, J. S., Mushotzky, R. F., & Arnaud, K. A. 1992, ApJ, 395, 275 
Dong, X. Y., & De Robertis, M. M. 2006, AJ, 131, 1236 

Elvis, M., Risaliti, G., Nicastro, F., Miller, J., & Puccetti, S. 2004, ApJ, 635, L25 
Ferrarese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 
Ghisellini, G., Haardt, F., & Matt, G. 1994, MNRAS, 267, 743 

Iyomoto, N., Makishima, K., Fukazawa, Y., Tashiro, M., & Ishisaki, Y. 1997, PASJ, 49, 425 
Kartje, J. P., Konigl, A., k Elitzur, M. 1999, ApJ, 513, 180 
Krolik, J. H. & Begelman, M. C. 1988, ApJ, 329, 702 
Magdziarz, P. & Zdziarski, A. A. 1995, MNRAS, 273, 837 



- 11 - 



Marconi, A. & Hunt, L. K. 2003, ApJ, 589, L21 

Matt, G., Guainazzi, M., & Maiolino, R. 2003, MNRAS, 342, 422 

Oliva, E., Origlia, L., Kotilainen, J. K., & Moorwood, A. F. M. 1995, A&A, 301, 55 

Raymond, J. C, & Smith, B. W. 1977, ApJS, 35, 419 

Risaliti, G., Maiolino, R., & Bassani, L. 2000, A&A, 356, 33 

Risaliti, G., Elvis, M., & Nicastro, F. 2002, ApJ, 571, 234 

Risaliti, G., Bianchi, S., Matt, G., Baldi, A., Elvis, M., Fabbiano, G., & Zezas, A. 2005, 
ApJ, 630, L129 

Risaliti, G., Elvis, M., Fabbiano, G., Baldi, A., & Zezas, A. 2005, ApJ, 623, L93 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Starling, R. L. C, Siemiginowska, A., Uttley, P., & Soria, R. 2004, MNRAS, 347, 67 
Urry, C. M. & Padovani, P. 1995, PASP, 107, 803 

Weisskopf, M. C, Brinkman, B., Canizares, C, Garmire, G., Murray, S., & Van Speybroeck, 
L. P. 2002, PASP, 114, 1 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 12 - 



Table 2: NGC 1365 - Spectral fits 
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a : Column density in units of 10 22 cm -2 , obtained fitting a model with a free photon index 
power law. b : Same, with the photon index frozen to the average value, < T >= 1.68. c : 
Normalization of the power law, in units of 10~ 3 Photons s _1 cm -2 keV -1 . In OBS 2 this 
refers to the intrinsic emission producing the observed reflected spectrum. d : Fixed value. 
d : ratio between the normalizations of the reflection and transmission components. 



